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ABSTRACT 

We investigate the formation and evolution of the pseudobulge in "Eris" , a high-resolution TV-body 
+ smoothed particle hydrodynamics (SPH) cosmological simulation that successfully reproduces a 
Milky Way-like massive late- type spiral in a cold dark matter (ACDM) universe. At the present 
epoch, Eris has a virial mass M Ylv ~8x 10 11 M , a photometric stellar mass M* = 3.2 x 10 10 M , a 
bulge-to-total ratio B/T = 0.26, and a weak nuclear bar. We find that the bulk of the pseudobulge 
forms quickly at high redshift via a combination of non-axisymmetric disk instabilities and tidal 
interactions or mergers both occurring on dynamical timescales, not through slow secular processes 
at lower redshift. Its subsequent evolution is not strictly secular either, and is closely intertwined 
with the evolution of the stellar bar. In fact, the structure that we recognize as a pseudobulge today 
evolves from a stellar bar that formed at high redshift, was destroyed by minor mergers at z ~ 3, 
reformed shortly after, and weakened again following a steady gas inflow at z < 1. The gradual 
dissolution of the bar ensues at z ~ 1 and continues until the present without increasing the stellar 
velocity dispersion in the inner regions. In this scenario the pseudobulge is not a separate component 
from the inner disk in terms of formation path, rather it is the first step in the inside-out formation 
of the baryonic disk, in agreement with the fact that pseudobulges of massive spiral galaxies have 
typically a dominant old stellar population. If our simulations do indeed reproduce the formation 
mechanisms of massive spirals, then the progenitors of late- type galaxies should have strong bars and 
small photometric pseudobulges at high redshift. 

Subject headings: galaxies: evolution - galaxies: bulges - method: numerical 



1. INTRODUCTION 

The formation of realistic galaxies in cosmological sim- 
ulations has been a challenge over the past decade. 
The massive spiral galaxies formed in these simulations 
have traditionally been too centrally concentrated, with 
small disks and large bulges, a problem asso ciated with 
the catastrophic loss of angular moment um ( | Navarro fc| 
|Steinmetz]|2000 Scannapieco et al!]|2Q12[) . Recently, ow- 
ing to the increase in numerical resolution and the im- 
provement of sub-grid recipes for star formation and feed- 
back, several groups have reported successful simulations 
in ACDM of disk systems over a r ange of masses, from 
field dwarfs dGovernato et al.l 



20101) 



to massive 



dPiontek & Steinmetz 2011; Guedes et al.||2011 



Okamoto 20T2j), to ellipticals (3o lansson 



spirals 



Brook 



|et al.||2012| ). The increasingly better resolution ol these 
types ol simulations is now allowing us to study the de- 
tailed mechanisms involved the formation of galaxies and 
their different components. In the following, we focus on 
the formation of the bright, central component of disk 
galaxies. 

Possible scenarios for the formation of galactic bulges 
include early build-up via mergers at high redshift, sec- 
ular evolution due to the presence of a bar slowly bring- 
ing gas to the center and turning into a sphero id via 
the buckling instability or resonant thicke ning (e.g. |Raha| 
et al.||1991|lDebattista et al.||2004[ |2005| ), and fragmen- 
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tation via gravitational instability producing giant gas 
clumps that s piral-in towards the center via dynamical 
friction (e.g. |Noguchi| |1998| jlmmeli et al.||2003| |Bour-| 
|naud et al | 2007^ These distinct mechanisms lead to 
the formation of galactic bulges of different characteris- 
tics. "Classical bulges" are thought to form via mergers 
flNaab fc Trujillo||2006| |Hopkins et aT1|2010| ), which can 
rapidly produce spheroids by violently destroying disks, 
and are characterized by a steep increase in density to- 
wards the galaxy center. These bulges are best fit by 
high Sersic index (n > 2) profiles, resembling the sur- 
face brightness distribution of elliptical galaxies. "Pseu- 
dobulges' Vare_observed instead to have disk-like density 
profil es dCarollo et aI1[T997l [l998l ICarollo fc Stiavelli 



1998) and kinematics (IKormendy 19931), low Sersic in- 



dex (n < 2) profiles, and are thought to originate from 



disk material via secular evo lution (see a review by Kor- 
mendy fc Kennicutt Jr|2004 [ ) induced by spiral structure 



or bars (Combes fc Sanders 1981; Pfenniger & Norman 
1 990 ; [Combes fc Elmegreen|1993|llJebattista et al.|2 004 ; 
Athanassoula 20051) . While in general pseudobulges have 
younger stellar populations than classical bulges, massive 



late- type sp irals tend to have old pseudobulges (|Carollo 



et al. 2007). 



Traditionally, detailed studies of secular evolution have 
been confined to idealized isolated galaxy models, which 
by construction neglect the effect of cosmological accre- 
tion of matter, external perturbations by infalling satel 
lites, and the fact that the disk itself can grow a nd be de 



stroyed several times as a result of mer gers (e.g. Robert- 
son et al. |2006 Governato et al. | |2009|). No table excep- 
tions are the work of Curir et al. (2006, 2007) who studied 
bar instabilities in cosmological simulations at relatively 



2 



low ma ss resoluti o n (a f actor of 17 lower than Eris), and 
that of Okamoto (2012) who does not focus however on 
the interplay between the stellar bar and the formation 
of the pseudobulge, as we do here. Studies of isolated 
galaxy models have shown that the proper modeling of 
shocks that lead to the accumulation of gas and the sub- 
sequent formation of stars at Linblad resonances is nec- 
essary for the formation of realistic spiral structure, and 
that the amount of gas inflow is strongly dependent on 
the gravitational torques exerted by stellar bars. There- 
fore resolving the bar is crucial to studying the formation 
of the bulge. High gas concentrations can lea d to the de- 
struction of the bar (Frie dli fc B enz 199 3) due to the 
increase of chaotic stellar orbits ( |Hasan &J Norman 19901 
Hasan et al. |1993), generate nuclear bars (Friedli fc Mar-| 
tinet ||1993| ), fuel active galactic nuc lei activity (|Phinney| 
1994), and dissolve bars into bulges ( |Norman et al.|1 996). 
The amount of gas in the disk also influences the pattern 
speed of the bar by increasing the precessing frequency 
Q — k/2 ( Friedli fc Benzf l993), due to a steady tra nsfer 
of angular momentum from the bar to disk stars (Sell- 



wood||198ip . Bars may also trigger density waves i n the 
dark matter distribution ( |Weinberg fc Katz||2005| ) that 
can enhance and displace the dar k matterann ihilation 
signal from the centers of galaxies (iKuhlen ~et al. 2012 



In this paper, we use the Eris simulation (|Guedes et al. 



2011 ) to study the origin and evolution of the bulge com- 



ponent of a Milky Way-like late-type system, and its re- 
lation to the assembly of the disk, with unprecedented 
detail. Eris appears to be the first cosmological hydro- 
dynamic simulation in which the galaxy structural prop- 
erties, the mass budget in the various components, and 
the scaling relations between mass and luminosity are all 
consistent with a host of observational constraints. Be- 
cause of its high resolution and remarkable agreement 
with the observations, Eris provides then an excellent 
laboratory for understanding the mechanisms that lead 
to the formation of a pseudobulge in a single massive 
spiral whose internal structure is well resolved. 

2. THE ERIS SIMULATION 

2.1. Simulation setup 

Eris follows the formation of a M vir = 7.9 x 10 11 M 
galaxy halo from z = 90 to the present epoch in a 
Wilkinson Microwave Anisotropy Probe 3-year cosmol- 
ogy. This "zoom-in" simulation was performed using 
the TV-body + s moothed particles hydr odynamics (SPH) 
code GASOLINE flWadsley et al ||2004| ). The version of 
the code used for Eris includes Compton cooling, atomic 
cooling, metal licity-dependent radiative cooling at low 
temperatures ( |Mashchenko et aT| | 2006|), a prescription 
for blastwave supernova feedback ( jStinson et al.| [2006), 
and the ionizing effect of a un iform U V background us- 



ing a modified version of the Haardt fc Madau] ( |1996 ) 
spectrum. The target halo was selected to have a quiet 
merger history with no major mergers (defined as mass 
ratio > 1/10) at redshift z < 3 (four minor mergers oc- 
cur at z ~ 3, and two more at z ~ 1 and z ~ 0.6). 
The high resolution region was resampled with 13 mil- 
lion dark matter particles and an equal number of gas 
particles, for a mass resolution of raDM = 9-8 x 10 4 M 
and raspH = 2 x 10 4 M . The gravitational softening 
length was fixed to 120 physical pc for all particle species 



from z = 9 to the present time, and evolved as 1/(1 + z) 
at earlier times. Star particles form in cold gas that 
reaches a density threshold of n$F = 5 atoms cm -3 at 
the rate dp*/dt = 0.1p gas /tdyn, where p* and p gas are 
the stellar and gas densities, and td yn is the local dy- 
namical time. They are created stochastically with an 
initial mass m* = 6 x 10 3 M distributed following a 



Kroupa et al. (1993) initial mass function. Supernova 
explosions deposit an energy of 8 x 10 50 ergs and metals 
into the nearest neighbor gas part icles, and the heated 
gas has its cooling shut off following |Stinson et al.| p006). 
The use of a high threshold for star formation has the ef- 
fect of increasing the efficiency of supernovae feedback 
through the injection of energy in localized high-density 
regions. This results in t he formation of an inhomoge- 
neous interstellar medium ( |Governato et al. 2010 Guedes 
|et al.pOll) ) and the preferential remova l of low angular 
momentum material (Br ook et al.||20lT ). 



2.2. A realistic late-type massive spiral 

At z = 0, Eris is a late- type spiral galaxy of virial mass 
M vir = 7.9 x 10 11 M and virial radius R vir = 239 kpc, re- 
solved with A/dm = 7xl0 6 ,iV gas = 3xl0 6 ,7V* =8.6xl0 6 
dark matter, gas, and star particles, respectively. To 
apply photometric analysis techniques to the simulation 
data, we use the radiative transfer code SUNRISE (|Jon-| 



sson 



2006] ) in post-processing to generate mock images 
of the galaxy in a range of filters. Eris' i-band abso- 
lute magnitude and total photometric stellar mass are 
Mi = -21.7 and M* = 3.2 x 10 10 M o , respectively. A 
two-dimensional decomposition of the SPSS i-ba nd im- 
age using the GALFIT package (|Peng et"aL 2002) shows 
that Eris has an extended disk of exponential scale length 
Rd = 2.5 kpc, a pseudobulge with S ersic index n = 1.4 , 
and a bulge-to-disk ratio B/D=0.35 ( Guedes et al.|20lT ). 
The present-day u — r = 1.52 mag color of Eris is con- 
sistent with the colors of g alaxies that host pseudob- 
ulges (Drory & Fisher 2007). Eris's surface brightness 
profile shows a down bending break at about five disk 
scale lengths, as observed in many nearby s piral galaxies 
flPohlen fc Trujillo||2QQ6] |Bakos et al.||2008| ). The break 
radius rt, r = 10.7 kpc, the extrapolated central surface 
brightness fi- m = 19.9 i-mag arcsec -2 , and the break sur- 
face brightness /ibr = 23.5 z-mag arcsec -2 are all in excel- 
lent a greement with the observations ( |Pohlen fc Trujillo 
[2006] ). 

Ens falls on the T ully-Fisher relation as measured by 
|Pizagno et al.|(|2007[), on the locus of the Ssfr-^ht plan e 



2008) 



occupied by nearby spiral galaxies ( |Bigiel et al 
and on the stellar m ass-halo mass relation at z = 
(Behroozi et al. 2010). Its rotation curve is consistent 
with observations of blue hori zontal-branch ha lo stars in 
the Sloan Digital Sky Survey ( |Xue et aT]|2008[ ). Its stel- 
lar disk is currently forming stars at a rate, 1.1 M yr _1 , 
comparable to the value inf erred for the Milky Way 
(Robitaille & Whitney 2010). Its cold gas reservoir, 
1.9 x 10 u M(?), is in good agreement with valu es inferred 
for the Milky Way ( |Nakanishi fc Sofue||2003| . Cold gas 
extends to six stellar disk scale lengths, and is pierced by 
kpc-size holes produced by SN explos ions as observed in 
nearby spirals ( Boomsma et al.||2008 ). 

Twin simulations of the Eris galaxy that include metal- 
dependent radiative cooling at high temperature (T > 
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Fig. 1. — Evolution of Eris from z = 4 (t oy) to z = (b ottom). Left panel: composite i-band, V, and far UV rest-frame images generated 
using the SUNRISE radiative transfer code ( Jonsson 2006). Right panel: Stellar surface mass density contours of the inner 5 comoving kpc 
of the galaxy, in face-on and edge-on projections through 10 kpc slices. The bar forms together with the proto-disk at z ~ 5, is destroyed 
by minor mergers at z ~ 3, reforms soon afterwards, and finally weakens from z < 1 to the present. 
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10 4 K) and exp licit diffusion of heat and metals among 
SPH particles (Wadsley et al. 2008) have been shown 
to reproduce quantitatively the proper ties of the circum- 
galac tic medium of galaxies at z ~ 3 dShen et al.| 2011 



2012]) , which reflect the complex exchange ol matter, en- 
and metals between galaxies and their surround- 



ergy 
ings. 



3. DEFINING THE PSEUDOBULGE 



Throughout this paper, the use of the term "pseudob- 
ulge" is dictated by the low Sersic index of the inne r 
regions of the galaxy ( |Kormendy fc Kennicutt Jr| 2004 ) . 

Before formally decomposing the galaxy in its mam 
components, it is instructive to take a qualitative look 
at the Eris galaxy and its inner regions across cosmic 
time. Figure [I] shows the photometric evolution of Eris 
(left panel) and its inner regions (right panel) in face- 
on and edge-on portraits. It is these inner regions that 
we identify as the "pseudobulge" , without trying to sep- 
arate out the bar-like component of the potential. On 
the left panels we show composite rest-frame z,V, and 
far UV imag es of Eris from z = 4 to z = generated 
by SUNRISE (Jonsson||2006|), and on the right we show 



iso-density contours of the stellar distribution within 5 
kpc of the center. The disk-like structure that character- 
izes the pseudobulge is present already at z = 4, and the 
stellar disk grows around this initial structure as time 
goes on. A stellar bar is present at this redshift, and is 
destroyed at z = 3, an epoch characterized by a series of 
minor mergers (see left panel). Between redshifts 2 and 1 
the bar reforms, but by z = the bar is reduced to a very 
small nuclear bar only a few softening lengths in size, and 
has left behind a flattened spheroidal pseudobulge (see 
Figure [T]) . 

In oraer to identify a subset of particles that are rep- 
resentative of the z = pseudobulge, we decompose the 
galaxy based on its photometric and kinematic proper- 
ties by selecting particles that (i) lie in the region where 
the light profile is dominated by the Sersic profile, and 
(ii) are not kinematically cold. 

Figure [2] (top) shows the i-band surface brightness pro- 
file of the galaxy at z = and the photome tric decom- 



posit ion obtained with the GALFIT software ( Peng et al. 
2002). The resulting parameters are: an integrated bulge 



magnitude Mi = —19.3, a B/D = 0.35, a Sersic index 
n = 1.4 projected face-on, an effective radius R e = 0.68 
kpc, and an exponential disk scale length R s = 2.5 kpc. 
The pseudobulge dominates the light profile within a ra- 
dius of ~ 2 kpc, or three effective radii. To isolate the 
warmer pseudobulge component and remove stars that 
are kinematically cold (and most likely belong to the 
thin disk), we decompose the star particl es kinemati- 
cally according to their angular momentum ( |Abadi et al.| 
2003; |Scannapieco et aT]|2009[ ) as follows. We define 
the quantity e = J z /J c , where J z = J Z (E) is the an- 
gular momentum of a star of total energy E along the 
z direction; e is calculated on a centered and aligned 
frame of reference, where the xy plane coincides with 
the plane of the galactic disk. The maximum angular 
momentum of a particle of total energy E is given by 
J c = V c * r, where V c = y^GM(< r)/r is the circular 
velocity. In this decomposition scheme, stars with low 
angular momentum occupy the region around e = and 
belong to the pseudobulge and stellar halo, while stars 
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Fig. 2. — Photometric and kinematic definition of Eris' pseudob- 
ulge at z = 0. Top panel: Eris' surface brightness profile in the i- 
band (black points), and its best-fit decomposition into Sersic bulge 
and exponential disk given by GALFIT. The dotted line shows the 
radius below which the pseudobulge component dominates the light 
profile. Middle panel: Angular momentum distribution of all stars 
(dotted line), of stars located within 3R e = 2 kpc (solid black line), 
and of the pseudobulge stars with R < 3R e and e = J z /J c < 0.8. 
Bottom panel: Surface density profile of the kinematically decom- 
posed pseudobulge (orange dotted line) and thin stellar disk (blue 
dotted line), as well as of all star particles within 20 kpc from the 
center. 



with e ~ 1 have nearly circular orbits and belong to the 
thin, centrifugally-supported disk. 

The resulting stellar angular momentum distribution 
is shown in the middle panel of Figure [2] The dotted 
line shows the distribution of all particles within 16 kpc 
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Fig. 3. — Origin of pseudobulge stars. Top panels: Stellar age distribution of the present-day pseudobulge (left) and thin disk (right). 
Star particles that have been accreted ("non in-situ star formation") are shown with the green curve. Middle panels: Formation radius 
(black solid line) and current radius (dotted blue line) of pseudobulge (left) and thin disk (right) stars. Bottom panels: Merger history of 
the pseudobulge (left) and thin disk (right). The y-axis shows the formation redshift, the x-axis the formation coordinate in the frame of 
the host galaxy. 



from the center of the galaxy, showing the prominence of 
the thin stellar disk at e ~ 1 and warmer components at 
lower values of e. It is clear that there is a continuum in 
the angular momentum distribution going from the disk 
to the pseudobulge and stellar halo (the peak at e = 0). 
This is consistent with the fact that the central pseu- 
dobulge is effectively the low angular momentum tail of 
the disk rather than a truly distinct component (see also 
Dubois et al. 2012). The solid line shows the angular 
momentum distribution of all particles within r < 2 kpc, 
and the dashed yellow line shows the same distribution 
after applying an angular momentum cut at e < 0.8, 
which excludes thin disk particles. Thus, in our defi- 
nition, the pseudobulge spans a wide range of angular 



momenta, having a hot component at e ~ but also a 
tepid, disk-like component at intermediate angular mo- 
menta. The surface density profile of the kinematically 
derived pseudobulge and thin disk are shown on the bot- 
tom panel of figure [2] 

In summary, the pseudobulge is defined as consisting 
of particles star particles that lie within a radius r < 2 
kpc and are kinematically warm with e < 0.8. Particles 
with e > 0.8 are assigned to the thin disk. 

4. FORMATION AND EVOLUTION OF THE 
PSEUDOBULGE 

Having defined a subset of particles that belong to the 
present-day pseudobulge, we proceed to study its forma- 
tion and evolution. 
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4.1. Origin 

The main mechanism for the assembly of Ens' pseu- 
dobulge is by no means secular and it is not distinct from 
the formation path of the inner disk. Figure [3] shows the 
distributions of age and formation radius, and the merger 
history of the present-day pseudobulge and thin disk. 
Pseudobulge stars form in an early burst at z ~ 3 — 5, 
while thin disk stars have a spread of ages. Accreted 
stars (green histogram) are identified as particles that 
were born in a subhalo other than the main progenitor. 
In order to track the origin of these particles, we use the 
AMIGA Ha lo Finder (AHF, |Gill et al.||2004| |Knollmann 



Redshift 
1 0.5 



fc Knebe]|2QQ9[ ) to identify the host galaxy and subhalos 
from z = 10 onwards. We find that only 4% the pseu- 
dobulge mass was accreted, and that the bulk of this 
component formed in situ. Similarly, 98% of the thin 
disk's mass formed in situ. 

The middle panel shows the distribution of formation 
radius inform and present-day radii R now of the pseudob- 
ulge (left panel) and the thin disk (right panel) that 
formed in-situ. Most of the pseudobulge stars formed 
within the inner 2 kpc of the center of the galaxy, and 
inward migration of stars in the region between 2 and 4 
kpc contributes only 13% of its final mass. The forma- 
tion of the thin disk occurs also mostly in-situ, although 
11% of all stars that formed within r = 4 kpc migrate 
outwards to 15 kpc, possibly via the resonant s cattering 
of stars by transient spiral arms (Ros kar et al.||2008| ). 

The formation histories of the pseudobulge and the 
thin disk (bottom panels) are remarkably similar, in the 
sense that mergers do not contribute to the formation of 
the pseudobulge more than they contribute to the forma- 
tion of the thin disk. In our simulations then, the pseu- 
dobulge is the building block of the rest of the galactic 
disk, and does not form a posteriori from instabilities in 
it. As a result, the age distribution in the plane of the 
disk is a continuum from the inner to the outer regions, 
as shown in Figure [5] 




Time [Gyr] 

Fig. 4. — Top panel: evolution of the amplitude of the m=2 
Fourier mode of the stellar density distribution. Middle panel: 
Evolution of the circular velocity and the velocity dispersion within 
2 comoving kpc from the center. Bottom panel: Evolution of the 
gas mass enclosed within different radii as a function of redshift. 
Large inflows followed by outflows occur during minor mergers at 
z > 4, at z = 3 and z = 1. Dotted gray lines represent minor 
merger events at z < 5. 



4.2. Evolution 

While in isolated galaxy models a massive stellar disk 
is in place in the initial conditions, in hierarchical galaxy 
formation the disk forms and grows over time, is at least 
partially destroyed during mergers, re-grows after them, 
and suffers continuing external perturbations by infalling 
matter, in particular satellites. All such phenomena play 
an important role in redistributing mass and angular mo- 
mentum in the disk, contributing to the formation of a 
central, more concentrated component that we can rec- 
ognize as the "bulge". At late times, when the galaxy 
is no longer harassed by mergers, less violent (although 
not strictly secular) processes begin to play an important 
role in the evolution of the galaxy. 

Secular evolution is not a single mechanism, but rather 
a collection of different phenomena that can in principle 
take place even on dynamical timescales as a result of 
rapidly changing conditions both internal and external 
to the disk, as expected in the hierarchical assembly of 
galaxies. These processes have been s tudied predomi- 
nantly with isolated galaxy simulations (Debattista et al. 
2006). The formation of a stellar bar is the first step in 
the process. The bar can then buckle into a kinemat- 
ically hot spheroid as a result of bending modes, and 



drive an inflow of gas, resulting in new star formation 
taking place place in the center. The strength of the bar 
and therefore the effectiveness of buckling and inflows 
depend on the relative proportions of gas and stars, with 
more gas tending to weaken the bar (Has arTet al. | |1993[) 
and p reventing buckling to first order ( [Deb attista et al 
2006). Both the reshaping of pre-existing stars by inter- 
nal instabilities and the new central star formation can 
contribute to the build-up of the bulge from disk mate- 
rial. The Eris simulation has a mass and force resolution 
comparable to some of the reference studies of bar insta- 
bilities and secular evolution published in t he literature 



(Debattista et al. 2006 Athanassoula 
added complexity of the cosmologica 



2008[ ) despite the 
It thus 



context. 

allows to address in detail the role of these different pro- 
cesses traditionally considered part of secular evolution. 

4.2.1. The Role of Minor Mergers 

Although only a small fraction of the mass of the pseu- 
dobulge is contributed by accreted satellites, minor merg- 
ers aid to its evolution by destroying the bar and redis- 
tributing angular momentum in the disk. 

Figure [I] shows that a star forming disk-like structure 
with a strong bar mode is already in place at z = 4 and 
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Fig. 5. — Map of the mass- weighted age distribution of the galaxy 
projected face-on through a 2 kpc slice. In terms of stellar age, the 
transition from the pseudobulge to the inner disk is smooth. Some 
key features are the spiral spiral structure contained within the 
pseudobulge and a weak, younger nuclear bar within it. 



evolves significantly afterwards. In order to quantify the 
strength of the bar as a function of redshift, we calculate 
the amplitude of the m Fourier mode of the stellar den- 
sity distribution in polar coordinates E(i2, 0,£) at each 
simulation output as follows: 



A m (R,t) = 



1 



f 

Jo 



E(R,e',t)e- im6 'de' 



The phase is given by 

4>(R,t) = tan" 1 



L^(A™)J 



(1) 



(2) 



where 5ft (A m ) and ^s(Am) are the real and imaginary 
components of A m , respectively. When a strong bar is 
present, the phase of the m = 2 remains constant over 
the radial extent of the bar i^bar, and changes drastically 
beyond it. Therefore the strength of the bar can 
be measured as the maximum amplitude within a radius 
Rbar- In cases where the bar is weak, we measure A2 and 
the maximum amplitude within a fixed radius R = 2 kpc. 
The top panel of Figure [4] shows the evolution of A 2 as 
a function of redshift. Several minor mergers (shown as 
dotted vertical lines) at z ~ 3 destroy the bar, which 
reforms quickly by z = 2 and survives for 1 Gyr. After 
the last minor merger at z = 1 the bar steadily weakens, 
becoming negligible at z = 0. 

It is important to stress that a bar forms as soon as 
the inner disk forms at high redshift (z >4), as demon- 
strated by the amplitude of the m = 2 mode at those 
early epochs. The disk therefore develops a concentrated 
inner disk (a "bulge") with a low Sersic index as typ- 
ical in present-day pseudobulges already from the very 
beginning and on fast, dynamical timescales. Such con- 
centrated inner disk, which coincides with the bar, un- 




FlG. 6. — Evolution of the e = J z /Jc distribution as a function of 
redshift for all stars within the optical radius. Disk stars populate 
the region around e = 1 (co-rotation) while stars that occupy the 
regions around e = form part of the spheroid, which includes 
pseudobulge and halo stars. Particles with intermediate angular 
momentum can be associated with either the pseudobulge or the 
thick disk according to their binding energies. The first structure to 
form has low angular momentum and dominates the distribution 
down to z ~ 3, at which point the disk component has grown 
massive enough to become dominant. 

dergoes continued evolution (Figures [I] and [I]). What we 
witness is the action of the processes normally included 
in secular evolution but happening in a dynamical fash- 
ion; indeed there are continued external tidal triggers, 
resulting in a sequence of repeated, fast developing non- 
axisymmetric instabilities interrupted by minor mergers 
rather than a slow growth of non-axisymmetric modes 
giving rise to a bulge-like component only after several 
billions of years. 

The evolution of the average velocity dispersion, a, of 
the central regions (Figure ffl is tightly interconnected 
to the merging history of the galaxy and the evolution 
of the bar. As the protogalaxy begins to form and the 
bar is violently destroyed by minor mergers, the velocity 
dispersion increases, reaching a peak of 200 km s _1 . After 
z = 3 the central velocity dispersion is only marginally 
affected by mergers and remains nearly constant from 
z = 0.6, to the present time. This is because late mergers 
are more strongly tidally disrupted by the disk and most 
of the remnant material does not become part of the 
pseudobulge, whereas at early times the bulge dominates 
the protogalaxy are is more vulnerable to mergers. 

Thus, the weakening of the bar at z < 1 is not accom- 
panied by an increase in velocity dispersion as would be 
expected if the buckling instability was the major driver 
of the formation of the pseudobulge. This is consistent 
with studies that show that the buckling instability re- 
quires long-lived strong bars coupled with high central 
stellar densities to ensue, whereas in Eris the bar weakens 
and becomes more compact at low redshift due to angu- 
lar momentum losses via dynamical friction and bar-dark 
matter halo interactions. 
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Fig. 7. — Evolution of the photometric properties of the galaxy 
as function of redshift. The top panels show the evolution of the 
photometric B/D, B/T and D/T ratios obtained by fitting a Sersic 
profile t o the bulge and an exponential profile to the disk with 
Galfit ( |Peng et al.|2 002). The bulge component includes the con- 
tribution ot the bar). The evolution is non-monotonic since z = 4. 
The second panel shows the evolution of the Sersic index n. At 
early times the bulge has a Gaussian-like profile (n ~ 0.5) which 
steepens at late times to reach a maximum value of only n ~ 1.5, 
suggesting that there is never a classical bulge in Eris. The bottom 
panel shows the growth of the scale length of Eris' disk. There is 
a remarkable similitude in the growth rate of the Sersic index and 
the scale length, suggesting that the formation of the disk and the 
pseudobulge are tightly interconnected. 
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Fig. 8. — Left: Surface mass density map of the gas within 5 kpc 
of the center of Eris at z = 0. Right: Origin of the inner gas disk; 
surface mass density distribution of the gas particles identified to 
be in the region r < 2 kpc at z = traced at z = 1. 95% of the 
gas that today is in the inner 5 kpc was at r > 50 kpc at z = 1. 



4.2.2. The Role of Gas Inflow 

In Eris, there is steady inflow of cold gas at z < 1, that 
originates from accreting cold flows entering the halo and 
feeding the galaxy, and continues within the galaxy itself 
down to its nucleus. It contributes to the assembly of the 
outer disk, to replenishing the gas reservoir at smaller 
radii, and to sustaining a rate of star formation of 1- 
2 M yr _1 during this time. Figure H (bottom panel) 
shows the total gas mass enclosed within 0.3 to 20 kpc. 
While the gas mass within the galactic disk only increases 
by a factor of ~ 2 from redshift one to zero, the steady 
inflow of gas contributes to the weakening of the stel- 
lar bar suggested by the decreasing power of the m = 2 
mode (top panel) by increasing the central baryonic con- 
centration thus destabilizing box orbits, and contributes 
to the build up of the pseudobulge with new star forma- 
tion but without increasing its stellar velocity dispersion. 
The fact that the inflow of gas has a cosmological origin 
is shown in Figure |8j where we have traced back the gas 
that is at r < 2 kpc today (right panel) to its location 
at z = 1 (left panel). Gas that is located today within 
2 kpc from the center, was funneled to the center via 
and dense cold gaseous filaments at high redshift. Mi- 
nor merger events, shown as dashed lines in Figure [4j 
with peaks in the gas mass distribution and followed by 
outflows, do not contribute significantly to the gas mass 
accretion. 



In addition to increasing the velocity dispersion, mi- 
nor mergers also redistribute angular momentum within 
the disk. Figure [6] shows the distribution of angular mo- 
mentum of all particles within the optical as a function 
of redshift. The first structure to form is composed of 
mostly low angular momentum material. By z ~ 3.5 the 
disk begins to dominate the distribution. The disk and 
spheroid grow together, except at z = 1 when a minor 
merger heats up the inner regions, destroys the bar, and 
redistributes angular momentum from high to interme- 
diate values. The inset shows a detail of this angular 
momentum redistribution before (green) and after (or- 
ange) the minor merger, showing that there is a transfer 
of angular momentum from the disk (e ~ 1) to the pseu- 
dobulge and thick disk (e < 0.8). 



4.2.3. Photometric Evolution of the Bulge and the Disk 

Figure [7] shows the evolution of the photometric prop- 
erties of the galaxy, obtained by running GALFIT on rest- 
frame i-band images produced by SUNRISE in a range 
of redshift s. Each image is fit with Sersic index and an 
exponential disk profile, and we do not attempt to de- 
compose the bar separately. The top panels show the 
evolution of the B/D, B/T, and D/T ratios. These ra- 
tios do not evolve monotonically; their value is dictated 
by the evolution of the disk, the bulge, and the bar. As 
shown in the kinematic decomposition (Figure |6|, the 
bulge component forms first and dominates the light dis- 
tribution until z ~ 3. The B/D ratio fluctuates in re- 
sponse to the strength of the bar, and at z < 1, the B/D 
ratio and the Sersic index n increase as the large-scale 
bar gradually dissolves leaving behind a fairly compact 
pseudobulge with a nuclear bar. Since the bar forms dy- 
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namically immediately after the disk arises, the B/D of 
the galaxy is within a factor of 2-3 from the final value 
since the beginning with only rare instances, such as af- 
ter the last major mergers, in which the bar/pseudobulge 
are temporarily destroyed. Instead of increasing the cen- 
tral cusp, mergers have the opposite effect here, at least 
at high redshift, where they destroy the stellar bar effec- 
tively lowering the photometric bulge-to-disk ratio. The 
scale length of the disk, doubles since z = 2 due to 
a steady gas inflow that builds the disk inside-out. Al- 
though the pseudobulge forms during an early burst of 
star formation, the Sersic index of the inner region of the 
galaxy is never high, implying that at birth the pseudob- 
ulge stars are kinematically tepid, not hot. 

5. SUMMARY AND DISCUSSION 

We have studied bulge formation in a late-type disk 
galaxy. The bulge that forms has all the characteristics 
of a pseudobulge, consistent with the fact that late-type 
galaxies, including the Milky Way, have pseudobulges. 
However, while it is normally believed that such bulges 
are the result of a slow growth from internal dynamical 
instabilities, termed "secular evolution" , we have shown 
that in our simulations the process is fast and dynami- 
cal. This is in very good agreement with observations of 
pseud obulges in Milky Wa y- like galaxies in the local uni- 
verse ( |Carollo et al.||2QQ7[ ). Bars develop on a dynamical 
timescale as soon as the disk begins to assemble. This, 
and not fragmentation into clumps, seems to be the main 
manifestation of disk self-gravity and internal dynamical 
evolution at the galaxy masses considered here. The re- 
cent discoveries of clumpy galaxies at high z, and the 
simulations developed to explain them, are compatible 
with this picture as long as one accepts that disk frag- 
mentation becomes a mode of internal evolution only at 
much larger disk mass, presumably corresponding to a 
halo mass in excess of 1O 12 M , as the Toomre parame- 
ter is easily broug ht below unity by accretion of gas an d 



dwarf galaxies ( Governato et al.||2010 ). 



disturbances (e.g. Agertz et al. 2009 Genel et ah]|2010|). 



The pseudobulge is part of the low angular momentum 
tail of the disk material rather than being a separate 
component, not surprisingly since it is disk material that 
was assembled at early times in the initial disk assembly 
phase (due to early instabilities) or that was brought 
to the center by gas inflows at later times. Therefore 
there is a age distribution of stars in the pseudobulge, 
which is not young but rather contains a substantial old 
component, contrary to common objections against the 
internal instability/secular evolution picture of bulge 
formation. If this scenario is correct, progenitors of 
late- type disks, which are the majority of present-day 
galaxies, are thus expected to have strong bars and 
small "photometric pseudobulges" (i.e. indicated by 
low Sersic index inner stellar profiles) already at z > 3. 
Such bars are a significant fraction of the galaxy size 
in length, more than in the present-day galaxies. At 
masses significantly lower than those considered here, 
supernovae feedback might maintain a larger gas fraction 
and prevent the formation of a massive, bar unstable 
stellar disk u ntil much later tim es, producing younger 
pseudobulges (Brook et al. 2011), or no bulge at all in 



Our results can be summarized as follows: 

• In our simulations, the bulk of the pseudobulge's 
mass forms early (z ~ 4), in situ, and rapidly (2 
Gyr) from non-axisymmetric instabilities such as 
bars, not secularly. While the disk undergoes sev- 
eral episodes of such instabilities it never fragments 
into massive clumps as proposed in other scenarios 
of early bulge formation. Minor mergers contribute 
only 4% of its final mass. 

• The formation path of the pseudobulge and of the 
thin disk are indistinguishable in terms of merger 
history, suggesting that the inner disk and the 
pseudobulge are not distinct components. And 
since in our simulations the age profile distribution 
is continuous from the inner to the outer regions, 
the pseudobulge is simply the building block of the 
galaxy and of the inside-out disk formation scheme. 

• We find that the evolution of the pseudobulge is 
tightly intertwined with the evolution of the bar. 
Present-day pseudobulge particles formed early in 
a bar-like configuration. The bar mode was de- 
stroyed at z ~ 3 by minor mergers, but reforms 
quickly by z ~ 2. The steady increase in gas mass 
due to the presence of the bar and the redistribu- 
tion of angular momentum within the stellar distri- 
bution at z ~ 1 lead to the gradual dissolution of 
the bar, leaving behind a compact flattened bulge 
and nuclear bar, and a steady increase of the Sersic 
index from n = 0.5 at z = 4 to n = 1.4 at the 
present epoch. 

Finally, we note that the pseudobulge formation pro- 
cess shown here is not peculiar to our simulation code 
and star formation prescr iption and parameters. Recent 
work by |Okamoto| ( |2012[ ) shows that pseudobulges that 
form in galaxies of the mass of Eris and higher form via 
similar mechanisms as shown here in simulations of only a 
factor of ~ 2 lower resolution, lower star formation den- 
sity threshold, and strong kinetic feedback giving rise 
to hydrodynamically decoupled winds. The early and 
rapid formation of galactic bulges in cosmological simu- 
lations suggests that attempts to produce massive bulge- 
less galaxies, s uch as those that seem t o dominate in the 

should focus on 



Local Volume ( Kormendy et al. 
suppressing this early burst of s 



2010]), 
:ar formation by intro- 



ducing, f or instance, a n additional mode of early stellar 
feedback |Brook et al.] fl2012| ) . 
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